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Abstract 

An exact solution for a high speed deterministic traffic model with open boundaries and 
synchronous update rule is presented. Because of the strong correlations in the model, 
the qualitative structure of the stationary state can be described for general values of 
the maximum speed. It is shown in the case of f max — 2 that a detailed analysis of 
this structure leads to an exact solution. Explicit expressions for the stationary state 
probabilities are given in terms of products of 24 x 24 matrices. From this solution an 
exact expression for the correlation length is derived. 
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1 Introduction 



One dimensional driven diffusive processes have proven to be an interesting playground for 
the study of non-equilibrium behaviour fl]-|5j . Of great interest is the fact that the study of 
many stationary state properties has come within reach of exact analytical methods since 
the solution of the asymmetric simple exclusion process (ASEP) with random sequential 
update and open boundaries |6|-[|. An important analytical tool in the study of these 
diffusive systems is the matrix product method, that appeared earlier in the study of 
lattice animals [ItJ and the ground states of antiferromagnets |TT|,[T2|. Its use in |J for 
the ASEP has boosted a lot of research on a variety of diffusion models, among which are 
for example the ASEP with other updates |TB|^TU|, multi-species models multi- 



lane traffic [g(J and the partially asymmetric exclusion process |31|-|34|. For recent reviews 
of many of the exact results for the ASEP see |35|-37|. 

It has been shown for different dynamical update rules, that the stationary state of a 
stochastic model can always be written as a matrix product [j38j-f40ll , although no proof 



is given for the synchronous update. This mere fact by no means solves the problem of 
finding the stationary state, but it provides a basis for a systematic study via the repre- 
sentation theory of non-linear algebras [^5|,|2B],[3^]. In almost all cases studied so far, the 
algebra has been quadratic, which is peculiar to systems with nearest neighbour interac- 
tions only. Only for the synchronous ASEP, non trivial representations of an algebra of 
higher degree have been used, i.e. either quartic |17j], where the matrices depend on one 
site, or cubic fl8| in the case where the matrices depend on two sites. 

Recently the asymmetric exclusion process with next nearest neighbour interactions 
has been studied by various methods |4l|] , but in the case of open boundaries exact results 
have been obtained only on a special line. In an initial attempt to find exact stationary 
states for models with long range interactions and open boundaries, a deterministic high 
speed asymmetric exclusion model is studied. Particles are allowed to hop over more 
than one lattice spacing per time step, and they enter the system at the left and leave at 
the right. Furthermore, the system will be subject to a synchronous dynamical update 
rule. For such dynamics the correlations are the strongest, which is not only interesting 
from a physical point of view, but actually helps solving the problem. The correlations 
are so strong that one can describe the stationary state qualitatively in a simple way by 
identifying so called Garden of Eden states p2[|43| . 

The exact stationary state is given in matrix product form for the case where particles 
may hop over two lattice sites. The matrices depend on three sites and it is shown that 
they should satisfy an algebra of at least quartic degree. The model is solved by making 
an Ansatz for the form of the matrices based on the qualitative observations for the 
stationary state. It is then shown that the submatrices in this Ansatz must satisfy certain 
relations for the matrix product state to solve the stationary master equation. A solution 
of these relations is found with the help of explicit calculations for small systems. 

This paper is organised as follows. The model is defined in detail in Section 0. The 
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exact stationary state is calculated in Section [5] and some results on the phase diagram 
are discussed in Section |j. 

2 Definition of the model 

In this paper we study a one dimensional asymmetric particle hopping model where 
particles in the bulk hop to the right. Particles may enter the system at the left and 
leave at the right. In the bulk all particles will move with their maximum possible speed, 
which is either given by the speed limit f max , or it is given by the distance to the next 
particle to avoid collisions. There will be no stochasticity in the bulk and particles always 
achieve their maximum possible speed instantaneously. In the case of periodic boundary 
conditions this model is known as the deterministic Fukui-Ishibashi model ||44| , for which 
some exact results are known |)45|| . 

With open boundaries, particles will be allowed to enter the system on the first t> max 
sites and may leave the system from the last v max sites. The choice of boundary conditions 
can have a profound influence on the behaviour of the system. If, for example, particles 
were allowed to enter only at the first site, the density profile for t> max > 1 in the free 
flow phase would show a strong sublattice dependence. Moreover, in this case the system 
would not be able to reach its maximum possible flow, since particles would have to wait 
an extra time step due to the synchronous update before the first site is unblocked. In 
the case of t> max = 2, the specific boundary conditions we will use here are similar to 
those of the random sequential model A of j|(J. If the two sites at the left boundary are 
empty, a particle can enter on the second site with probability a 2 , and on the first site 
with probability — a 2 ). The sites remain empty with probability (1 — q;i)(1 — a 2 ). 
If a particle is already present on the second site, a probability a 3 is given for a particle 
entering on the first site. At the right boundary, a particle at the last site will leave 
the system with probability If the last site is empty, but a particle is present on 
the penultimate site, it will leave the system with probability (3 2 . in terms of Boolean 
variables 7j that have the value 1 for a particle and for a hole, the dynamical rule for 
the bulk can be written as, 

where the prime denotes time incremented by one and a = 1 — r. At the boundaries 
additional Boolean variables aand/3 are used that have time averages equal to a and (3. 
At the left boundary the rules are, 

r{ = tit 2 + di(l - a 2 )<ri<T2 + a30-iT2, (2) 
T ' 2 = ri a 2 T 3 + r 2 r 3 + a 2 a x a 2 , (3) 

and at the right boundary we have, 

T' L = T L _ 2 L _ X L + (1 - fc)T L + (1 - hYL-X°L. (4) 
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The currents are defined by the continuity equation, 

r- =Ti+ji-i -ji, (5) 

and are given by, 

jo = a 3 criT 2 + (1 - (1 - «i)(l - a 2 ))a 1 a 2 
ji = na 2 + a 2 a x a 2 

ji = + Tj_i(7jCrj + i (6) 

JL = Ar L + P 2 r L -icr L . 

For technical convenience we will put (/3i) = (/3 2 ) = /3, and (di) = (a 2 ) = (63) = ot in the 
rest of this paper. All arguments however hold for the more general case. The calculation 
becomes more cumbersome and one has to discriminate between even and odd sublattices. 



3 The stationary state 

In the following discussion, the relative weight of a particular configuration {t±, . . . ,tl} 
in the stationary state will be denoted by P(t± : . . . ,t l ). Once all relative weights are 
determined, the normalisation Zl can be calculated via, 

Z L = £P(7i,...,T L ). (7) 

M 

To derive some general conclusions about stationary state, it is helpful to consider the 
extreme cases a — 1 and j3 — 1 first. 

3.1 Free flow 

Because we are considering a deterministic model, spontaneous jams do not occur. Jams 
will only build up from obstacles at the right boundary. Pure free flow configurations 
are obtained by removing these obstacles, i.e., (3=1. The dynamical rule at the right 
boundary then becomes, 

t' l = r L _ 2 o L _ x o L . (8) 
Since there will be no jams in the stationary state, its bulk dynamics is given by, 

< = Tm, (9) 

and it follows that the master equation can be written as, 

P(ti, . . . ,t l ) = ^2pf(riT 2 T 3 )p { (r 2 T 3 T A )P(T 3 , . . . ,r L , //,//), (10) 
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where, 

p f (000) = (1 - a) p f (100) = a 

Pf(010) = 1 Pf(OOl) = 1 (11) 

Pf^rr'r") = otherwise. 
This equation can be solved by the Ansatz, 

L-2 

P(ri, . . . ,t l ) = i?(r L _!r L ) JJ p f (riT m r i+2 ), (12) 

and we find, 

P(00) = 1, R(10) = R(01) = a. (13) 

3.2 Jammed flow 

Pure jammed flow configuration are obtained by setting a — 1. From the dynamical rules 
it follows that in this case configurations with the sequence 000 in it do not occur in the 
stationary state. This means that the bulk and left boundary dynamics may be replaced 
by the simple rule, 

< = r m . (14) 

The master equation for this case is, 

P(n, . . . , t l ) = ^pj(r L _ 2 rL-ir L )P(/i, n, . . . , r L -i), (15) 
n 

where, 



(16) 



Pj(100) 


= p 




Pj(010) 


= P 




Pj(001) 


= 1 




Pj(H0) 


= P 




Pj(101) 


= 1 - 


p 


Pj(0H) 


= 1 - 


-p 


Pj(lll) 


= 1 - 


p 


Pj(rr'r") 


= 


otherwise 



Again this equation can be solved by a simple Ansatz, 



L-2 



P(ti, ...,t l )= L(tiT 2 ) JJpj(rir i+1 r i+2 ). (17) 



i=i 



In this case we find, 

L(00)=/5 2 , L(10) = L(01) = P, L(ll) = l-p. (18) 
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3.3 The general case 

As mentioned already above, spontaneous jams will not occur since we are considering a 
deterministic model. Following a similar line of reasoning as in [43 , this can be deduced 
from the microscopic dynamics (|T])-(||). 

i The sequence 1100 can only arise from the same sequence shifted by one lattice unit, 

(-TjTi+l 0-1+2^+3)' = T-T i+ iT i+2 <Ji + 3<7i + 4- (19) 

Since, 

(^-3^-2^-1^)'" ~ {r L _ 2 T L _ l( r L )" ~ {t l . x t l )' = 0, (20) 

it follows that configurations with the sequence 1100 do not occur in the stationary 
state. 

ii Similarly, a sequence 10100 can only arise from the same sequence shifted by one 
lattice unit or from a sequence with 1100 in it, 

{ T iO~i+l T i+2°~ i+?jO~ i+i)' — ( r i-2°"i-l + ^-l "*) ( r i+l cr i+2 T i+3 ~i+4°~i+5) 

+ T-a i+ iTi + 2Ti +3 <Ji + 4<Ji +5 . (21) 

Since, 

{tL-4<7L-3TL-2<7L-1<7l)'" ~ (TL-20L-1TL)' = 0, (22) 

it follows with the previous observation that also configurations with the sequence 
10100 in it do not occur in the stationary state. 

It thus follows that each configuration can be divided into three parts. The first part 
is a free flow part where there are at least two holes between successive particles. This 
part ends at site / which denotes the last site of the last 000 sequence of a configuration. 
The dynamics for this part is given by, 

t[ = r,_ 2 , 3 < i < f. (23) 

The third part starts at site j which denotes the first site of the first jammed configuration, 
i.e., a 11 or a 101 sequence, whichever comes first. This part is a jammed flow part where 
there are at most two holes between successive particles. For this part the dynamics is, 

r/ = r m , j<i<L-l. (24) 

In between these two parts there may be a sequence of 100's of arbitrary length. A general 
configuration may thus be written as 

71...T/(100) B T i ...71. (25) 
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This analysis can be performed in a similar way for models with higher f max - The free 
flow part will end with v max + 1 zeros and the intermediate part will consist of a sequence 
of blocks, where each block starts with a 1 followed by v max zeros. The jammed flow 
part starts with any of the local jammed configurations. These are those sequences where 
there are less than v max zeros in between two l's. 

The master equation for the stationary state can be written explicitly in this notation. 
In the case where the jammed flow starts with all pair it is given by, 

P(Ti . . . T/(100) n llr i+ 2 ...T L )= Pf ( T2 T3 ) Pf ( T2 T3 T4 ) p j ( T~L — 2 7~L — 1 7~L ) X 

P(r 3 ...T f (100) n+1 llr j+2 ...T L - 1 ) 



+ J2P(n... ^(100)^001(100)^11^+2 . . . Tjj- 



p=0 

n 



+ ^ P(r 3 . . . r / (100) p 010(100)"- p llr, +2 . . . r L _ x ) 

p=0 



(26) 



A slightly different equation is obtained when the jammed flow starts with a 101 sequence, 

P(ri . . . T/(l00) n 101r :?+ 3 • • • T L ) = Pf(TiT2T 3 )p{(T2T 3 T4 : )p i (TL--2TL-lTL) X 

P(T3...r / (100r +1 10lT J+3 ...T L _ 1 ) 



+ P ( r 3 • • • T/(100) p 001(100) n - p 101 



T j+3 ...TL-l, 



p=0 

n 



+ ^ P(r 3 . . . r / (100) p 010(100) n -P101r 3 - +3 . . . t l ^) 
p=o 

+P(r 3 . . . T f {W0) n+1 01lT j+3 . . . r L _i)] . 



Similar equations are obtained when / and/or j are close to the boundary. 

To solve (p6|) and (p7 ) we will employ the powerful matrix product method 
relative probabilities for the stationary are written as, 



L-2 



P(t 1 ,...,t l ) = (C(t 1 t 2 )\1[M( )\n{r L ^T L )), 



(27) 
The 

(28) 



i=l 



Because of the specific form of each configuration, the following triangular form for the 
matrices M. suggests itself, similar to the t> max = 1 case fl8fl, 



M(tt't" 



Pi (rr'r")F S(tt't") 
Pi (rr'r")J 



(29) 
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where Pi{tt't") and p^tt't") are defined by ( |TT| ) and (|T6| ) respectively. The matrices F 
and J are yet to be determined. While for f max = 1 they are just scalars given by F = (3 
and J = a, they will be more complicated in the present case. The matrices S{tt't") will 
solve the dynamical equations for the bulk. They are defined on the interface only and 
5(000) = 5(110) = 5(101) = 5(011) = 5(111) = 0. A similar decomposition as in (||) 
will be used for the boundary vectors, 

(C(t 1 t 2 )\ = ((C f (t 1 t 2 )\,(C,(t 1 t 2 )\), (30) 

and likewise for \TZ(tl-iTl)) ■ 

To make the following more transparent, we will use the notation 5i = 5(100), 52 = 
5(010) and 5 3 = 5(001). Upon substitution one quickly concludes that ( PE| ) and QZ7| ) are 
equivalent if, 

S 2 J = aF(S 2 + (1 -a)S 3 ). (31) 

Let us assume that this relation is satisfied and concentrate on (|26|) . Substituting (|29|) in 
( |26D one finds that, 

(£ P (00)| = (£ P (10)| = (£ F (01)| = (£ P |, 

|ftj(00)> = l^j(lO)) = |ftj(01)> = l^j(ll)) = \Kj), (32) 
and that the bulk matrices must satisfy, 

n-l 

a p [3 2 ( n - p -VF 3p+2 ({3 2 S 3 J 2 + f3FS 2 J + F 2 S 1 ) J 3 ^)" 1 + a n F 3n+2 S 3 J = 

p=0 

n 

a p f3 2{n - p) F 3p (f3 2 S 3 J 2 + f3FS 2 J + F 2 S l ) J 3 («-p) +1 + a n+1 F 3(n+1) 5 3 

p=0 

+ (i _ a )2 (1 a Pp^-p)F 3p+2 S 3 J 3(n - p)+1 

p=0 
n 

+ (1 - a)(l - (3) ^a p [3 2{n - p) F 3p+l (/35 3 J + F5 2 ) J^-pH 1 . (33) 

p=0 

The requirement that the four sums in (^) cancel term by term leads to the following 
equation, 

F 2 ((3 2 S 3 J 2 + (3FS 2 J + F 2 5 x ) = (3 2 ((3 2 S 3 J 2 + f3FS 2 J + F 2 5 x ) J 2 

+ (3 2 {l - a)(l - P)F ((1 - a)F5 3 + (3S 3 J + F5 2 ) J 2 , . (34) 

For the the remaining terms in (|33"D to cancel, the following equation must be satisfied, 

F 2 S 3 J = aF 3 S 3 + ((3 2 S 3 J 2 + (3FS 2 J + F 2 S l ) J 

+ (1 - a)(l - p)F ((1 - a)FS 3 + f3S 3 J + FS 2 ) J. (35) 
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Altogether we get three relations, ([34]) and (pop. These can be rewritten as, 

S 2 J = aF(S 2 + (l-a)S 3 ), 

= (3{a(3S 3 + S 2 )J + FS 1 (36) 
aF 2 S 3 J = a(3 2 S 3 J 3 + F((l-a)(l-/3)S 2 + a(3(l-a-(3)S 3 )J 2 
+ a 2 F 3 S 3 . 

Besides these bulk relations, there are boundary relations that follow from considerations 
of cases where / and/or j is close to the boundary. They are not particularly illuminating 
and are listed in appendix It is important to note that if we manage to solve these 



boundary relations such that (|32|) is also satisfied, a solution of (p6|) then ensures station- 
arity of the matrix product state (|28| ) for arbitrary system sizes. This also means that 
fl32D and (|36|) enable us to extrapolate knowledge of small systems to arbitrary large ones, 



which helps us to find a solution of the relations we have obtained. To find a representa- 
tion for the matrices F and J, we employ the usual strategy for these type of problems: 
to consider explicit solutions for small systems and to try to find relations between them. 
Using the Ansatz fl29|) for the particular form of the matrices, we then deduce that the 
following algebraic relations hold, 

= F 3 - (3(1 - a - a(3)F 2 - a(3 2 (2 - a)F - a 2 (3 4 

= (F - (3)(F + a(3)(F + a(3 2 ) -a(3 2 (l -a)(l - 0)F. (37) 

while J satisfies, 

= J 3 - a (2-a)J 2 + a 2 (3(l-a-a(3)J + a A (3 2 

= (J -a)(J -a(3)(J -a 2 (3) + a(l-a)(l- (3) J 2 . (38) 

Considered as polynomials, each of these equations has three solutions. These may be 
thought of as being the eigenvalues of F and J respectively. One thus finds a three dimen- 
sional representation for F and J for which we have to check that it is compatible with 
the other relations. This is indeed the case and two examples of explicit representations 
for all objects are given in appendix |B[ For a,j3 ^ 0, F will have non-zero eigenvalues 
and thus is invertible. It is then found that ( |38D is satisfied by choosing J = —a 2 j3 2 F~ 1 . 

Although the solutions of the cubic equation ( p7|) are awkward expressions in terms 
of a and (3, they are all real for < a,j3 < 1. Let A n denote the eigenvalues of F and 
H n = —a 2 (3 2 /X n the eigenvalues of J. They can be written in the following way, 

\ n = a + pcosl J, (39) 

fd>'- (2n + 1)tt\ . . 

fin = 6 + p'sin r K — '—). (40) 



with Ai < A 2 < < A 3 and fi 3 < < fi\ < fi 2 for < a, (3 < 1. Here, a, b, p, p', and 4>' 
are real functions of a and (3 defined in appendix |B.l . 
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4 Results 



In this section some exact results concerning the phase diagram are calculated. In the 
case of the purely free flow and jammed flow phases, the current and density profiles can 
be calculated easily. For general values of the boundary rates the correlation length is 
calculated and it is shown that it diverges on a special line in the phase diagram. 



4.1 Free and jammed flow 

We have seen that the stationary state in both the free flow case (f3 = 1) and the jammed 
flow case (a — 1) is a simple product. This means that correlations are absent and the 
density profile is flat in both cases. The values of the current and the density are easily 
calculated and given by 

a . 2a 

pF = iw JF = rw (41) 

pj = j> = -y^t- (42) 

It follows that in the free flow phase j-p = 2p F . As expected all particles move with their 
maximum speed f max = 2. In the jammed flow phase, the fundamental diagram is given 
by J J = 1 — Pi- I n this case all holes move with their maximum speed which is equal to 
1. Note that the values of the two currents are equal if 2a = (3(1 + (3). We will see in the 
next section that this line is the coexistence line in the general phase diagram. 

4.2 The general case 

To facilitate summing over the stationary state probabilities, the local free and jammed 
flow configuration weights pf(tt't") and pf(tt't") are collected into matrices Pf and P } 
with elements, 

(P flJ ) Tr , TV , = Rj(rrV'), (43) 
= otherwise. 

Likewise, the matrix S and vectors (L\ and \R) are defined, 

S TT ,, TlT „ = S(tt>t"), (L\ tt , = (C(tt% \R) tt , = \R{tt')). (44) 

The advantage of this notation is that for example the normalisation Z L can be written 
compactly as, 

Z L = J2 P(n, ...,t l ) = (L\M L ~ 2 \R), (45) 

{r} 
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where M is a 24 x 24 matrix given by, 



To perform the calculations, it is worth mentioning the following two intermediate results, 

((P^jriRj))* = -^-(1,1,1,1)® (^+ 2 , ^ +2 , Ms +2 ) , (48) 
where the discriminant A is defined by, 



A = (A 1 -A 2 )(A 2 -A 3 )(A 3 -A 1 ). 



(49) 



After some laborious manipulations, the normalisation is then found to be, 



(3\l-a)(l-(3) 
a(2a- (3(1 + (3))/\ 



a: 



73(1 - a)(l + 2a) ]T ^±^(\ l+1 - A,+ 2 )Af 



i=l 



+ (i _ ^ £ - „ +2 )ft 1+1 

^— ' (/ij - a) 2 



i=l 



(50) 



where A3+j = Aj and /ij+3 = //j. This expression is well defined for all a and (3. To see 
this, the only non-trivial case we have to consider is 2a = (3(1 + (3). All other explicit 
poles in Q50D are equivalent to (1 — a)(l — (3) = 0. At 2a = (3(1 + /3) the solutions 
Aj simplify dramatically. In particular Ai = — (3 2 while A 2 and A3 are the roots of a 
quadratic equation. Furthermore, we find that A 3 = /i 2 and A 2 = /i 3 . These values imply 
that the expression between brackets in ([50]) has a zero that precisely cancels the pole at 
2a = (3(1 + (3). 

Unfortunately the calculations with the present notations are still rather intricate and 
laborious and so far no other explicit expressions for the general case have been obtained. 
The phase behaviour of the model however is similar to that of the case v max = 1 P^|18| , [4T 



From the explicit form of the normalisation it is clear that the correlation length will be 
determined by ratio of its largest contributions and we find, 



£ F = In — 



in 



for 2a<p(l + p), 
for 2a > (3(1 + (3). 



(51) 
(52) 



There is a low density phase for 2a < (3(1 + (3) and a high density phase for 2a > (3(1 + (3). 
In each case, the bulk density will have the free and jammed flow value respectively. At 
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the boundaries there will be exponential corrections of which the correlation lengths are 
given by fl5lD and fl52|). The curve 2a = [3{1 + (3) is a coexisting curve on which the 
correlation length diverges. The instantaneous density profile is a shock profile resulting 
in an average linear profile. Across this curve, the average bulk density has a jump of size 
Pj ~ Pf- 

As we have seen above, the locus of the coexisting curve is obtained by equating the 
values of the current for the two extreme cases: the free and jammed flow phases. This 



seems to be a general feature of ASEP's and supports the ideas of Kolomeisky et al. [46 



for the case of discrete time and parallel update. In the present model however, the 
dependence of the correlation length on a and (3 does not decouple, as is the case for 



1 [|18j and the random sequential ASEP ||. It is therefore quite amazing that the 



locus of the coexisting curve still can be obtained by a simple mean field analysis. 



5 Conclusion 

An exact solution for the stationary state of a deterministic traffic model with v max = 2 is 
presented. Apart from the absence of symmetry due to the lack of a particle-hole duality, 
the phase diagram is qualitatively similar to that of the case with t> max = 1, as expected. 
This solution might be a first step towards an exact solution of a realistic traffic model. 

The stationary state is presented in a matrix product form, where the matrices depend 
on three sites and are 24 dimensional. The matrix product method has been shown to 
work extremely well in those cases where the matrices are either infinite dimensional or 
of small finite dimension. Its shortcomings are obvious when the matrices are of large 
finite dimension and the eigenvalues become solutions of polynomials of high degree. As 
in the present case, it will still be a tedious technical exercise to derive exact expressions 
for expectation values and correlation functions from the exact solution. 

Although many relations between matrix elements have been given, no proper matrix 
algebra has been derived. It will be interesting to find this underlying algebra, which at 
least should be of degree 4 as suggested by equation (|56]). This algebra may provide more 
convenient ways of deriving expectation values and correlation functions than the method 
used in this paper. The appearance of cubic roots however will remain. 

An obvious and interesting extension of the model will be to include stochasticity in 



the bulk hopping rates. This can be done as in the Fukui-Ishibashi model |)44] , but more 
interesting perhaps will be the aggressive driver Nagel-Schreckenberg model | 4"?|] , which 
is closer to real traffic but might be still simple enough to be analytically tractable. 

As a final remark one should mention that because of the long range interaction it is 
unlikely that these models are integrable in the sense that there would be and underlying 
Yang-Baxter relation. The paradigmatic ASEP however is integrable since it is closely 
related to the integrable XXZ spin chain. It would be interesting to compare the matrix 
product ground states of non-integrable systems with those of integrable ones |H7],[|H . 
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A Boundary relations 

In the cases where / < 3 extra boundary relations are needed for (|28| ) to be the stationary 
state. The following equations can be deduced from the master equation in those cases, 

• / = 3. 



a(C F \FS 3 J = a/3 2 (£j(01)|J 3 + (l-a)(l- (3)(C F \S 2 J 2 

+ a(3(l -a- (3)(C F \S 3 J 2 + a 2 (C F \F 2 S 3 . (53) 



/ = 2. 



a(3{\ - /3)(£j(01)| J = /3(£j(10)| J + {C ¥ \S±. (54) 



a(C F \S 3 J = af3 2 (£ j(01)|J 3 + a(l - a) (1 - (3) ((C F \S 2 + (1 - a) (C F \S 3 ), J 

+ a(3{\ - /3)(Cj(01)\J 2 + a 2 (C F \FS 3 . (55) 



/ = !■ 



/? 2 (£j(01)| J 2 + (£ F | {(3S 2 J + FS X ) = (3 2 {l - (3) ((1 - a)(£j(01)| + (£j(10)|) J 2 , 

(56) 

(£j(01)| J = a(C F \S 3 + (1 - p) ((1 - a)(£j(01) + (£j(10)|) J. (57) 



/ = o. 



(A(H)| = ^(l-/5)(£j(01)|, (£f(H)| = 0, (5J 



/5 2 (£j(11)| J 2 = /3(£j(10)| J 2 + (£ F |5iJ. (59) 
Similarly, when j > L — 1 extra relations are needed. These are, 
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S 2 \Kj) = a(l - a)\K F (00)). 



J = L 



F 2 S 3 \TZj) = aF 3 \1Z F (00)) + {(3 2 S 3 J 2 + (3FS 2 J + F 2 S 1 )\1Z J ) 

+ (1 - a)(l - /3)F((1 - a)FS 3 + (3S 3 J + FS 2 )\K 3 ). 

• j = L + 1 

(1 - (1 - a) 2 )F 2 |ft F (00)> = (3((3S 3 J + FS 2 + (1 - a)FS 3 )\n 3 ). 

B Representations 

B.l The diagonal representation 

In this subsection an explicit representation is given in which F and J are diagonal. 
Ai, A 2 and A 3 denote the three solutions of the equation, 

A 3 - 3aA 2 - 36A - a(3 2 c = 0, 

where, 

a = (Ai + A 2 + A 3 )/3 = /3(l-a-a/3)/3. 

b = -(AiA 2 + A 2 A 3 + A 3 Ai)/(3/3 2 ) = a(2-a)/3. 

c = A!A 2 A 3 //5 2 = a 2 /? 2 . 

Then, F and J are given by, 

F = diag{Ai, A 2 , A 3 }, J = -a^diagjAr 1 , A,' 1 , A 3 

The matrices <Sj are given by, 

/ -Ai - A 2 Ai + A 3 

51 = a 2 /9 2 Ai + A 2 -A 2 - A 3 

\-Ai - A 3 A 2 + A 3 

/ -AiA 2 AxA 3 

5 2 = (1 - a) I AiA 2 -A 2 A 3 

\— AiA 3 A 2 A 3 

AiA 2 + a/5 2 -A!A 3 -a/5 2N 

-AiA 2 - a{3 2 A 2 A 3 + a(3 2 

AiA 3 + a/3 2 -A 2 A 3 - a{3 2 
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The corresponding boundary vectors are given by. 
(£ F | = a(3 2 (l ~ a)(l - (3) 



2/i N/1 a\ I ^1 ^2 A3 



Ai-/3' \ 2 -p A 3 -/3, 
(£f(11)| = 0, 

(£j(10)| = ((Ai - /3)(A 2 - A 3 )((2a - l)Ai + a 2 /?), (A 2 - /?)(A 3 - A x )((2a - 1)A 2 + a 2 /3), 
(A 3 - /3>(Ax - A 2 )((2a - 1)A 3 + a 2 /?)) (67) 
(£j(01)| = a/3 2 (l-a)(l-/3)^ Al(A2 ~ A3) A2 ^ 3 " Al ) ^(Ax - A 2 ) 

(£j(H)l = ^(l-/3)(£j(01). 



Ax + a/3 2 ' A 2 + a/3 2 ' A 3 + «/3 2 



and the right boundary vectors are given by, 

l*j> = -^(Ar 2 ) A 2 2 ,A 3 - 2 ), 

I^f(OO)) = -^(Ax(A 2 -A 3 ),A 2 (A 3 -Ai) ) A 3 (Ai-A 2 )) ) (6£ 
|ft F (10)> = |fl F (01)) = |i2 F (ll)> = 0, 



where the discriminant A is defined in (|49|). 

The solutions of ( |63|) are all real for < a, /3 < 1 and can be written in the following 
form, 

f(b + 27m\ . . 

A n = a + pcosl , (69) 

where, 

p = 2 v / a 2 + /3 2 &, <P = arctan (A/C) . (70) 

The discriminant A and C are given by 



A = v / 27p 6 /16 - C 2 , C = 3v / 3(2a 3 + /? 2 (c + 3a6)). (71) 
Similarly, the solutions \i n = —a 2 (3 2 / \ n of ([38|) can be written as, 

,„^ + ,-sin( ^-< 2 3 " + ^ ), (72) 

where, 



p' = 2Vb 2 - a 2 a, 0' = arctan (C'/A) . (73) 
C = 3V3(3 2 {2b 3 -a 2 {c + 3ab))/a 2 . (74) 
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B.2 A simple representation 



Although the eigenvalues of F and J, or equivalently the roots of (|37j ) and (p8[), are 
awkward expressions, an example of a representation with simple matrix elements is given 

by, 

/ 1 a \ /0 a/3 

F = P 1-a -a/3 a(l - a) , J = a /3 1-a 1 - /3 | . (75) 

\ 1-/3 -a/ \0 1-a 1 

The corresponding representations for Si, 1S2 and S3 are given by, 

5! = --^f-= U/3(l-/3) a(l-a)(l-/3) a(l - /3) 2 

1 ^ W(l-a) a(l-a) 2 a(l-a)(l-/3) 



/3 2 



1-/9 



1-a 1-a 0\ /0 1 - a 1 

a(l-a)(l-/3) , S 3 = (3 -a/3 



(76) 



a(l — a) a(l — a) 2 0/ \0 a(l — a) a. 

The left boundary vectors in this representation are represented by, 

(£ F | = (l,l-a,0), 
(£j(10)| = ^(/3(l-a),(l-a)(a + /3-a/3),a(l-/3)) 

(£j(00) = 0, (£j(01)| = (0,l-a,l) (77) 

(A(n)l = 1(1-^X^(01), 

and the right boundary vectors are given by, 
\K 3 ) = (0,1-/3,1), 

n (78) 

\n F (io)) = |^f(01)) = |^f(11)) = 0, |^ F (00)) = -(l,0,a). 

a 

This representation is convenient for calculations for small system sizes. For larger system 
sizes it is more useful to use a representation in which F and J are diagonal. The price 
one has to pay is that the matrix elements will be more complicated because they are 
cubic roots. 
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